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Abstract 

We investigate the phenomenology of flavored dark matter (DM). DM stability is guaranteed 
by an accidental symmetry, a subgroup of the standard model (SM) flavor group that is not 
broken by the SM Yukawa interactions. We consider an explicit realization where the quark part 
of the SM flavor group is fully gauged. If the dominant interactions between DM and visible sector 
are through flavor gauge bosons, as we show for Dirac fermion flavored DM, then the DM mass is 
bounded between roughly 0.5 TeV and 5TeV if the DM multiplet mass is split only radiatively. In 
general, however, no such relation exists. We demonstrate this using scalar flavored DM where the 
main interaction with the SM is through the Higgs portal. For both cases we derive constraints 
from flavor, cosmology, direct and indirect DM detection, and collider searches. 
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I. INTRODUCTION 


The stability of dark matter (DM) is commonly assumed to be due to an exact discrete 
symmetry, Z n . This can either be imposed by hand or have a dynamical origin. Examples 
include R-parity in the MSSM [1], and flavor symmetries in the leptonic sector [2-5]. Here, 
we explore the intriguing possibility raised in Refs. [6, 7] that the stability of DM is due to the 
Z 3 center group of the global U| M = SU{3)q x SU(3)u x SU(3)d quark flavor symmetry. 
While (/f, M is broken by the SM Yukawa interactions, its subgroup Z 3 remains unbroken 
in the SM. More generally, it remains exact also in the presence of New Physics (NP), if 
the flavor breaking is of Minimally Flavor Violating (MFV) type, i.e. only due to the SM 
Yukawas. The lightest neutral state that is odd under Z 3 is therefore stable and is a DM 
candidate. This is the idea behind MFV dark matter [6-8]. 

Requiring MFV is sufficient, but not necessary. In this paper we formulate a general 
condition for flavored DM using flavor triality (see Eq. (3) below). For example, any spurion 
in the bifundamental of C/f. M leaves the above Z 3 unbroken. The flavor breaking can thus 
be quite far from MFV and still have stability of DM guaranteed by the flavor dynamics. 
To illustrate this point we consider the model of Ref. [9] where the flavor-breaking spurions 
have the form Yflj and are thus canonically not of the MFV type. In this model the SM 
quark flavor symmetry C/fU is fully gauged giving rise to flavor-gauge bosons, the Yukawas 
are promoted to physical scalar fields (flavons) transforming under flavor, and in addition 
there is a set of chiral fermions that cancel the anomalies in the flavor-gauge sector. 

Using this renormalizable model we show below that a thermal relic DM can be in a 
nontrivial representation of Qfl' 1 . There are two conflicting constraints on this setup. On 
the one hand, Flavor Changing Neutral Current (FCNC) constraints impose lower bounds 
on the masses of states in nontrivial flavor representations. On the other hand, a DM relic 
density consistent with observations requires large enough DM annihilation cross section so 
that some of these same particles need to be sufficiently light. Both of these requirements 
are satisfied for (P(TeV) DM mass. This is low enough that it may be tested by direct and 
indirect DM detection experiments and searched for at high-energy particle colliders. 

While the phenomenology of flavored DM models can be found in Refs. [7, 10-20], the 
construction of an explicit renormalizable model with inclusion of flavor-gauge interactions 
is new. Within our framework, the constraints on DM are more severe compared to a generic 
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Effective Field Theory (EFT) analysis [6, 8]. In particular, the flavor constraints from new 
fermionic states, and the fact that the vacuum expectation values (vevs) of the flavon fields 
need to reproduce the quark masses, makes the structure of the theory much more rigid and 
predictive. 

The paper is structured as follows. In Sec. Ill we derive general conditions for DM to be 
stabilized by the exact accidental flavor symmetry of the SM (flavor triality). An explicit 
realization of this possibility is introduced in Sec. Ill in the form of a model with fully 
gauged C?f. M . The resulting DM, flavor, and collider phenomenology is analyzed in detail in 
Sec. IV. We summarize our conclusions in Sec. VI, while more technical details of some of 
our computations are relegated to the Appendices. 

II. STABILITY OF FLAVORED DARK MATTER 

We start by formulating the general conditions required for flavored DM to be stable due 
to flavor triality. The SM exhibits a large global flavor symmetry U(3)q x U(3)u xf/(3)cX 
U(3)l x U(3)e in the limit of vanishing Yukawa interactions. In this paper we focus on the 
quark sector. This has the global symmetry x U{l)y x U(1)b x [7(1)pq. The three 
17(1) factors are the hypercharge, baryon number (5), and the Peccei-Quinn symmetry, 
respectively. The remaining semisimple group is £7f- M = SU(3)q x SU(3)u x SU(3)r>. The 
SM quarks transform under it as 


Ql ~ (3, 1 , 1), Ur ~ (1, 3,1), Dr ~ (1,1, 3). (1) 

The C/fN global symmetry is broken by the SM Yukawa terms 

Cy = QlHuuUr + Q L Hy d D R + h.c., (2) 

where H = ia 2 i7*. Cy is formally invariant under Qp if y u ^ are promoted to spurions that 
transform as (3,3,1) and (3,1,3) [21, 24, 71-75]. NP is of MFV type if y u . c i are the only 
flavor-breaking spurions also in the NP sector. 

The SM Yukawa couplings in Eq. (2) break [7(1)pq and break <7f, M to its center group 
, under which all three generations of quarks transform as 
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{Ur, Dr,Ql} —» U 2jt ^{Ur, Dr,Ql} ■ Ill the SM the Z^ DC2 is identical to a subgroup of 
V(1)b- This is no longer true in the presence of NP. In MFV for instance, Z^ D< ^ remains 
exact, while U{1)b can be broken, e.g., by dimension-9 operators [22] (see also [7]). 

The Z^ D< ^ may be the underlying reason for the stability of DM. To make this explicit it 
is useful to introduce Z£, a diagonal subgroup of Z^ DQ x Z\. Here, Z% is the center group 
of color SU(3) C , under which {Ur, Dr,Ql} —> e — i27r / 3 {f/ R , Dr, Q l } . All the SM fields are 
thus Z?{ singlets. In MFV NP Z£ is exact, so that the lightest Z£ odd particle is stable and 
can be a DM candidate [7]. 

We generalize this observation beyond MFV. To this end, we introduce the notion of 
flavor triality [7]. Consider a field X in the representation A" ~ (nW ttiq) x (rip, mfl ) x 
(rip, nip), where nf, rnf are the Dynkin coefficients of the corresponding SU(3)i group. We 
call flavor triality the quantity 

(rix — mx) mod 3, (3) 

where n x = rig + np + rip and mx = rrig + rn§ + nip. 

The basic requirements for flavored DM to be stable due to Z£ are the following. First 
of all, needs to be a good symmetry in the UV. Secondly, C/f, M needs to be broken only 
by spurions (<I>) with zero flavor triality, (n($) — m($))mod3 = 0. This ensures that Z\ is 
unbroken. (The spurions (<E>) need to be color singlets in order not to break color.) The 
lightest Z% odd state is then stable. If it is a color singlet it is a potential DM candidate. 
This also means that DM is in a nontrivial flavor representation with nonzero flavor triality, 
(n x — m x ) mod 3 fl 2 0. 

The above shows that models with flavored DM can deviate significantly from MFV. 
In particular, Z£ is not broken by a vev of any field that is in an adjoint or in a bifun¬ 
damental of Specifically, any function f(y u ,yd) leaves Z£ unbroken. More generally, 
additional flavor-breaking sources that transform as (8,1,1), (1,3,3),... may be present 
without spoiling DM stability. While the flavor structure of such NP models is not of MFV 
type in general, the stability of DM is still a consequence of an unbroken flavor subgroup. 
DM is in a nontrivial representation of the flavor group, leading to distinct phenomenology 
depending on the nature of the flavor breaking and on which flavor multiplet y belongs to. 

An important starting point in the above discussion was the assumption that is a 
good symmetry in the UV. This is most easily achieved, if t/fA is gauged. We explore this 
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possibility in the remainder of the paper. 


III. GAUGED FLAVOR INTERACTIONS AND DARK MATTER 


We gauge the full SM quark-flavor symmetry Q^. The fermionic sector is extended to 
cancel the anomalies of the new gauge sector. We use the model of Ref. [9] that allows for 
C?(TeV) flavor gauge bosons (FGBs). The SM Yukawas, y Uj d, arise from the vevs of new 
scalar fields transforming as 

Y u ~(3,3,l), Y d ~ (3,1, 3), (4) 

under The minimal set of chiral fermions that ensures anomaly cancellation of the new 
gauged sector is 

^uR ~ (3,1,1), *dR~( 3,1,1), (1,3,1), tf dL ~(l,l,3), (5) 


where the index L and R represents their chirality. Together with the SM fermions they, 
therefore, form vector-like representations of The SM gauge quantum numbers of 

^uR,^dR^uL,^dL are the same as for U R ,D R ,U R ,D R , respectively, i.e., they are SU(2) L 
singlets but charged under U(1)y- Because the new fermions are vector-like under the 
SM, e.g, T uR transforms like T uR under the SM, all anomalies in the SM sector cancel. 
Remarkably, with the above fermionic content also all mixed gauge anomalies cancel. In 
fact, one could also gauge two additional flavor diagonal f/(l)’s, U(l)u and U(1)d, as well 
as U(1)b-l, a possibility that we do not pursue further, but is discussed in Ref. [9]. 

The Yukawa and relevant mass terms in the Lagrangian are [9] 


^mass D XuQlH'&uh T X u ^ uR Y u ^/ uR T Yf u \I ! u lU R 

+ XdQ lH^ dR + X! d ^ dlXd^ dR + Md^dLD}I + ll.C., 


( 6 ) 


where \ u d are flavor-universal coupling constants and M U}d flavor-universal mass parameters. 
The mass terms in Eq. (6) mix the states T u l, u r and Ul,r forming mass eigenstates Ui and 
where i — 1,2,3 is the generation index (and similarly for down-quark states). After 
electroweak symmetry breaking the masses for the two mass-eigenstate sets are, in the limit 


m uf . > m Ui , [9, 23] 


v \ u Ad u 

ViKlYX 


m u[ ~ X u (Y u )i. 


(7) 
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The mass matrix for the FGBs, A\, A = Q,U, D and a = 1,..., 8 , is governed by the vevs 
of the Y u d scalar fields and the gauge couplings, [9, 23] 

(. M 2 AB ) ab = 1 -9 a9b Tr [(Y u ){X a , X h }(Y u ) ] ](5 AB 8 AQ -28 AQ 8 BU + Q O U)+U,u<* D,d, ( 8 ) 

with A a the Gell-Mann SU(3) matrices. The mass matrix is 24 x 24 dimensional. We 
denote the mass-ordered eigenstates by A m , m = 1 ,... ,24, where A 24 is the lightest one. 
The lightest gauge boson is found to be along the three diagonal A§ directions. This is a 
consequence of (Y u ) and (Y d ) being almost aligned and with very hierarchical eigenvalues, 
where the (Y u ) 33 and (Yd )33 are the smallest eigenvalues. 

The SM Yukawas, y u ,d, are generated after Y u j obtain a vev and the T, fields are inte¬ 
grated out. To first order in M Ujd /(Y Utd ), this gives 

A U M U XdM d 

y “ ~ kY) ’ s '~wi' 

The SM Yukawas, y u ^ d , are non-analytic functions of the spurions (Y Ujd ), which signals that 
the theory is strictly speaking not MFV. Analogously, the NP states, u ), d) and A m , have 
masses that are non-analytic in terms of the SM Yukawas. However, the low-energy ob¬ 
servables, with only the SM fields on the external legs can be MFV-like. If the M Ujd / (Y Ujd ) 
suppressed terms are kept in Eq. (9), the y U)d become more involved functions of {Y u ,d)~ l ■ 
These are analytic in (Y u ,d)~ l since the effects of NP states decouple in the (Y u j) —* 00 
limit. Similarly, the NP contributions to the low-energy observables C\ take the form 
SCi = F((Y U ) _1 , (Y rf ) _1 ) = F(y u ,y d ), with F,F analytic functions. One can thus expand 
5Ci = aiy u y\ + CL 2 (y u yi) 2 + b 1 y u yly d yl + • • •, where we assumed for illustration that the 
transition is due to the left-handed quark current. As long as there are no large flavor- 
conserving ratios, i.e., as long as (X U M U /X' u )/(X d M d /X' d ) l/|Wb|, the Taylor expansion 

can be truncated after a few terms (see Ref. [24] for a more detailed discussion). I 11 this 
limit, the low-energy effects are of the MFV type, suppressing the FCNCs to acceptable lev¬ 
els already for NP states at the electroweak scale. In a numerical analysis that we perform 
in Appendix A, we find that the expansion of the effective weak Hamiltonian in terms of 
SM Yukawa couplings can indeed still be performed for FGB contributions. 

Since (Y Ujd ) are in the bifundamental representation of the remains unbroken. 
As argued above the can be used to make flavored DM candidates stable. We consider 
two examples: i) a fermionic DM in a vector-like representation of that thermalizes 


7 




A QiA U: A D 


X 

Figure 1: Radiative corrections due to FGBs, Aq, A u , Ap, split the DM multiplet x- 

with the visible sector through FGBs, and ii) a scalar flavored DM, that interacts with the 
visible sector by exchanging FGBs as well as the Higgs. 


A. Fermionic flavored dark matter 


The DM in the first model is a massive Dirac fermion, y, in a vector-like representation 
of SU (3){j, 

(1,3,1), (10) 

so that no gauge anomalies are induced. Its mass term is 


£mass = m ° x XX- 


( 11 ) 


Since x is charged under Z£, the lightest member of the x triplet is stable. Note that we 
could also gauge a larger global group Qp M x SU( 3) x , with x transforming under SU(3) X . 
That we identify SU(3) X with SU(3)u is a model-building choice. 

The DM flavor triplet, y, is split by radiative corrections due to the exchanges of FGBs, 
see Fig. 1. In the -C tha limit, the DM mass splitting at one-loop is given by 


a 3 gfj o 

A m-y = - 

x 4 16tt 2 x 


- -Tr: 
3 


( 12 ) 


where A m x is a 3 x 3 matrix and so is 5 = A a (log M.\/^ 2 ) ab \ b . The FGB mass matrix 
M.\ is given in Eq. (8), while the a,b indices run only over the eight SU(3)u generators. 
The unphysical /j dependence cancels in the r.h.s. of Eq. (12). The Xu ' l = 1,2,3, mass 
eigenstates are obtained by diagonalizing the mass matrix A m x . And S is a function of 
Y,}Y U and YjY d Y d f Y u vev combinations, making the x mass eigen-basis slightly misaligned 
with respect to the up-quark one. The xi mass receives contributions from the heaviest 
FGBs (cf. section V). The lightest state is thus Xi , i- e -, with the predominantly up-quark 
flavor, while the heaviest is the top-flavored state, y 3 . 





8u 

Figure 2: Typical radiative splitting of the fermionic DM multiplet with A?B 3 i (Am 32 ) shown 
in red (blue) as a function of gjj, while all other parameters are kept fixed at gQ = 0.4, go = 0.5, 
M u = 600 GeV, M d = 400 GeV, \ u = 1, \' v = 0.5, \ d = 0.25, \' d = 0.3. 


In the numerics we use the exact one-loop expressions for the DM mass splitting, 

s = W m '° +8 A”. 


(13) 


Summation over FGB mass eigenstate indices m = 1,..., 24 and over a, c = 1,..., 8 is 
understood. The 24 x 24 dimensional matrix W diagonalizes the gauge-boson mass matrix 
and Bq i are the Veltman-Passarino functions. Typical values of the mass splitting as a 
function of gu are shown in Fig. 2. Denoting Amy = m Xi — m Xj , we see that Am 32 Am 3 i, 
so that the lightest state xi is split away from y 2 and X 3 3 with the latter approximately 
degenerate, m X2 ~ m X3 . This is very different from MFV DM [ 6 - 8 ] where the DM mass 
splitting is assumed to be expandable in the SM Yukawa couplings. In that case one has an 
approximate U( 2) symmetry for the first two generations giving m Xl ~ m X21 while the top- 
flavored DM, X 3 , is split away from the first two generations, and can be either significantly 
heavier or lighter. 

The relation Am 32 Am^i signifies that our flavored DM is non-MFV. The flavor gauge 
group SU(3)u is broken by the FGB vevs (Y u ). This breaking is larger in the first two 
generations. Since the quark masses are inversely proportional to (Y u ), this leads to an 
approximate global U(2)u symmetry in the quark sector. Such an approximate symmetry is 
not found for the radiative corrections to DM masses, m Xi . The DM multiplet has a chiral 
symmetry in the limit m x —>- 0 ensuring that the radiative corrections are proportional to 
m x and only log-dependent on FGB masses. The splitting does not vanish in the ( Y u ) — > 00 
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limit (or, equivalently, y u —* 0 limit), since in this limit the SU(3)u gauge group is still 
completely broken. Numerically, for m x ~TeV the splittings are Am .31 ~ (9(10 GeV) and 
A m 32 rv-/ (9(1 GeV) for gu ~ 0.4 and can be less than a pion mass for an order of magnitude 
smaller gu- 

The DM multiplet can be split more significantly if there is flavor breaking beyond ( Y u ), 
(Yd). As an example we consider an additional scalar held in the adjoint of SU(3)u, 4*c/ ~ 
(1, 8 ,1). The DM mass Lagrangian now reads 

£mass = ™\xlXr + A x Xl$uXr + h.c., (14) 


and yields DM masses that are split already at tree level, A m x = A x ($u)- We assume that 
($[/) is aligned with Y^Yjj. Then the two diagonalize in the same basis giving (9(1) splitting 
between all three members of the multiplet. The alignment is not needed in general, but 
does simplify our analysis. For the same reason, we also take the first state to be the lightest 
one, m Xl < m X2 3 . 

The Xi interact with the SM through FGBs. This also induces the decay of the heavier 
two states in the DM multiplet, y 2 , 3 , to Xi- We parametrize the relevant interactions with 


3 ($”) nXiYXiA 1 ; 

+(( si) U:m P R +&) u , m p L ) u,a- (is) 

+d,r((e d R ) k , m PR + (g d L ) il:m p,)d l A^, 

where Pr.l = ^(1 ± y 5 ) and k, l — 1 ,..., 6 . The couplings of Xi to the gauge bosons are 

(%h = (-lgvU'\VW m U) ii , (16) 


where U diagonalizes the m x mass matrix, Wm x U = rh x , and W diagonalizes the gauge- 
boson mass matrix. The explicit form of FGB couplings to exotic and SM quarks, (G'I'jr) kl m , 
can be found in Appendix A.2 of Ref. [23]. 

The partial decay width for Xi —• y XjQiQk is, neglecting hadronization effects, 


Q /\ 7T7,^ • ri 

r (^ ^ xmk) = [| 


9 

mi 


■(d) 


kl,m 


+ L R 


(17) 


where the sum runs over the FGB mass eigenstates m — 1,... ,24. Expression (17) is valid 
in the Ar/iy <g; rn Xi limit, and neglecting the quark masses. The above approximations are 
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valid for all the values of parameter for which the correct relic abundance is obtained and 
the FCNC, collider and direct DM detection constraints are satisfied. 

If the mass splitting is less than the pion mass the decay Xi —> XjQiQk is kinematically 
not allowed. The heavier Xi states then decay radiatively to \i Xjll■ For our purposes 
an order of magnitude estimate of the decay width suffices. The naive dimensional analysis 
estimate gives 

r (x. Xjii) ~ (iD 

(18) 



where Q u = 2/3 and Qd = —1/3 are the electromagnetic charges of up and down quarks. 
The sum over m runs over the FGB mass eigenstates, while the sum over / is over the SM 
quarks and exotic states, of mass m/ (for up, down and strange quarks this needs to be 
replaced with Aqcd)- 


B. Scalar flavored dark matter 

The second model has scalar DM, </>, in a fundamental representation of SU(3)u 

(1,3,1). (19) 

The main difference with the fermionic flavored DM from the previous subsection is that the 
scalar DM interacts with the visible sector via two different types of interactions. The first 
is its couplings to the FGBs, which is similar to the case of the fermionic DM. The second 
is a direct coupling to the Higgs 

£« = ( 20 ) 

For a thermal relic the DM annihilations proceed predominantly through the Higgs portal. 
The interactions via FGBs are subdominant except if ~ m a A /2 for some A a . Unless this 
is the case, the fact that the DM carries a flavor quantum number is exhibited only through 
the multiplicity of the states. 

After electroweak symmetry breaking, the DM-Higgs interactions are given by 

D A H (vh + u 2 /2 )0V , (21) 
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and the DM mass term is shifted by the Higgs condensate to give the DM mass of 

= ml +v 2 /2 . ( 22 ) 

The vevs of the flavons, (Y u ) and (Yd), split the DM multiplet at tree level through 

C D ki (<j>'\ a <j>)Tr(YZ\ a Y u ) + k 2 (^{ A°, A 6 }0)Tr(Yj{A o , A fe }Y n ). (23) 

The spectrum is also split by radiative corrections due to FGBs. These are quadratically 
divergent and proportional to the square of the FGB mass. In principle, it is possible to fine 
tune the tree-level and loop contributions to give almost degenerate DM flavor multiplet. 
However, given the hierarchical FGB masses, it is more likely that the DM multiplet is split 
completely, and only the lightest state is relevant for DM phenomenology. Depending on the 
signs of Ki in Eq. (23) the lightest 0 component can be either top-quark or up-quark flavored. 
We choose the latter option in the numerics for easier comparison with the fermionic case. 


IV. DARK MATTER AND NEW PHYSICS PHENOMENOLOGY 

We turn next to the phenomenology of the flavored DM models. We perform a scan over 
the parameters of the models and show that the lowest DM states, both for the fermionic 
DM, X :> and the scalar DM, 0, can be thermal relics. To make the scan numerically tractable 
we rely on several approximations in calculating the relic density, which we explain below. 
We also discuss the predictions for direct DM detection, and the constraints from FCNCs 
and collider searches. 

A. Scan results 

In the scan we fix A„ = 1 and vary A d G [1/(47t), 1]. The range is chosen with the expec¬ 
tation that one will be able to accommodate both the SM top and bottom quark Yukawas as 
well as satisfy electroweak precision constraints and direct t' and b' searches [9] . In addition 
we vary conservatively \' ud G [1 / (47t) 2 , 47t] , gQ t u,D £ [1 /(47t) 2 , 47t] , and M U)d G [0.2,20]TeV. 
To a good approximation, the variation of M u effectively compensates the fact that we do 
not vary X u as seen from Eq. (9). We have verified that further extending these parameter 
ranges does not extend the viable DM-model parameter space. For instance, the upper 
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Figure 3: The results of the scan for fermionic DM with radiative mass splitting (upper left 
panel), in the large mass splitting limit (upper right panel) and scalar (lower panel) flavored DM. 
Constraints from perturbativity (grey), t' (dark magenta) and dijet resonance (orange) searches, 
flavor bounds (light magenta), early-time cosmology (blue) and direct DM detection (brown) are 
consecutively applied. Allowed parameter points are denoted by green. For scalar flavor DM (right) 
we show the LUX and perturbativity bounds as two grey bands. The four benchmark points for 
fermionic flavored DM are denoted by a diamond, a triangle, a hexagram and a pentagram. 

ranges of gQ,u : D and \' u d already lie in the non-perturbative regime. To ensure perturbative 
control we require that all the FGB decay widths satisfy < 0.5 and that the ra¬ 
diative mass splitting for the fermionic DM is < 0.5 m x . This imposes upper bounds 

on gQ,u,D that are typically close to a/Ht. Similar constraints on X' u d are expected to follow 
from analogous considerations in the flavored Higgs sector, i.e., by requiring the total widths 
of the flavored scalars Y u ^ to be small compared to their masses. 

The results of the scan are shown in Figs. 3, 4, 5. Fig. 3 (upper panels) show the re¬ 
sults of the scan for fermionic DM model with radiative (left) and large tree-level mass 
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splitting (right). Fig. 3 (lower panel) instead shows the results of the scan for scalar DM. 
All the points shown in Fig. 3 give the correct relic DM abundance, Ddm- Different colors 
denote consecutively applied constraints. The grey points fail the perturbativity require¬ 
ment, T A m < 0.5 | Amy | < 0.5 m x discussed above. The points in brown are excluded 

by direct DM detection, the points in dark magenta by t' direct searches and the points 
in orange by dijet resonance searches. The flavor bounds exclude points in light magenta, 
while cosmological considerations - mainly from big bang nucleosynthesis - exclude points 
in dark blue. The green points are allowed by all constraints. In Figs. 4 and 5 we also show 
the points where it is not possible to obtain the correct relic abundance (denoted by light 
blue), and denote by dark red the points excluded by the combined direct-detection, collider 
and flavor constraints. 

For fermionic DM the observed relic abundance requires resonantly enhanced annihilation 
through s-channcl exchange of the lightest FGB, A 24 , see Fig. 6 (left). This leads to the 
correlation m Xi ~ 24/2 shown in Fig. 3 (upper panels). It is possible to obtain the correct 

relic abundance also if the DM mass is only approximately half of the lightest FGB (points 
away from the diagonal in Fig. 3 (upper panels)). These points require at least some of the 
couplings to be large and are excluded by flavor, collider, direct detection, or perturbativity 
constraints. For the allowed points the relation m Xi ~ m^/2 is satisfied to within a few 
decay widths of A 24 , i.e. to within (9(10%). The scalar flavored DM, on the other hand, 
predominantly annihilates through the Higgs portal, see Fig. 6 (right). There is thus no 
relation between and as seen in Fig. 3 (lower panel). 

In the remainder of this section we discuss how the various constraints on the DM model 
were obtained. 

B. Thermal relic 

For the calculation of the DM relic density we follow Refs. [25, 26]. To speed-up the 
numerical scan we work in the non-relativistic limit, using the freeze-out approximation. 
The details of the calculation are given in Appendix B. Among viable parameter points we 
choose four benchmarks that satisfy all other experimental constraints. For the benchmark 
points we verify the DM relic abundance calculation using the MadDM [27] package. We 
computed the required Feynman rules using the Feynrules [28] package. 
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Figure 4: The ratio of masses of the next-to-lightest to the lightest FGBs, for ra- 

diatively split DM multiplet (upper left panel), and for the large mass splitting limit (upper right 
panel), as functions of the DM mass, m xi , for the fermionic flavored DM. Lower panel shows the 
relative radiative mass splitting in the DM multiplet. The constraints due to perturbativity (grey), 
too large relic abundance (light blue), early cosmology (dark blue), flavor and direct bounds (dark 
red), are applied consecutively, leaving allowed points (green). 

1. Fermionic dark matter 


In the fermionic DM model the DM annihilation to quarks is dominated by s-channel 
exchange of the lightest FGB, A 24 , see Fig. 6 (left panel). The XiXi kijUj annihilation 
cross section is given by 


°(XiXi UjUj) ~ 


(g?)l s 1/2 (»+ 2'<) (W);, 24 + 


47 T 




Am 2 

Xt 


(s - m 2 A24 ) 2 + m\ 2A r 2 A24 ’ 


where 


SiU = (&±&)/ 2 . 


(24) 


(25) 
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Figure 5: The maximal decay time of the two heavy states in the DM multiplet as functions of 
DM mass (left) and the minimal mass splitting in the DM multiplet (right) for radiatively split 
fermionic flavored DM. The color coding is as in Fig. 4. 




Figure 6: The Feynman diagrams for the dominant processes in the DM annihilation for fermionic 
(left) and scalar (right) flavored DM. For scalar DM only one representative diagram is shown; other 
relevant final states include 66, cc, tt and tt,hh,ZZ (when kinematically allowed). 


a/s is the center of mass energy and T^24, are the decay width and mass of the lightest 
FGB, respectively. In Eq. (24) we have neglected quark masses; the full expression is given 
in Eq. (B8). The XiXi ~^ djdj annihilation cross section follows from Eq. (24) by replacing 
u —y d. The full decay width of the lightest FGB is the sum of all partial decay widths for 
kinematically allowed channels, 


r/4 


24 


UjUj) ~ 


m 


24 


47T 


(Sv) 2 a 24 + (dO 


2 

33 ,24 


(26) 


In the above expression we have neglected the quark masses for simplicity, with the full 
expression given in Eq. (B9). The rates for A 24 —> XiXi , djdj are obtained by trivial coupling 
replacements and by correcting the color multiplicity factors. 
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The correct relic abundance requires resonant annihilation, m x ~ 771,424/2, see Fig. 3 
(upper panels). This implies an upper bound on the DM mass through the following ar¬ 
gument. The thermally averaged DM annihilation cross section scales in the narrow width 
approximation as 


(av) oc 


9%4 


+ 0^,424/777,424) 


( 27 ) 


m A 2iT A 24 ' \ (y)2, 4 

Here, we used the approximate scaling for the FGB masses and decay widths, 777,424 ~ 
(Y)A 24 g A 2 i, T A 24 ~ (g A 24 .) 2 m A 24 ,. The ( Y) A 2i and g A 24 are, respectively, the projections of 
the Y Uj d vevs and gQ,u,D couplings onto the lightest FGB, A 24 . The DM relic abundance 
is Ddm oc 1 /(crv) oc (Y )^ 24 and thus depends predominantly only on the flavon vevs. Not 
exceeding the relic abundance puts an upper bound (Y) A 2 4 < 0 (few 100 GeV), almost 
independent of the DM mass. Since 777,424 ~ ( Y) A 24 g A 24 , and g A 24 < for the theory 

to be perturbative, this also sets an upper bound on the lightest FGB mass. This in turn 
implies an upper bound on the DM mass through the relation m x ~ 777,424/2. 

In the limit where only Xi contributes to the DM relic abundance we find, using the scan, 
an upper bound m xi < 10 TeV. The approximation is valid if y 2 ,3 states decay well before 
Xi freezes out (i.e. r X23 < 10 -11 s for m x ~ 1 TeV). For purely radiative DM mass splitting 
this is never the case (c.f. Fig. 5 ). Instead, if y 2 ,3 decay after decoupling, they will also 
contribute to the final DM relic abundance and one needs to sum all three contributions. In 
this case, the constraints on the mass spectrum become much more severe. In particular, 
in order for all y components to annihilate efficiently their masses need to be within a few 
decay widths away from the lightest FGB (LFGB) mass. This in turn implies that the 
(radiative) DM mass splitting has to be of the order of the LFGB width. Since the splitting 
increases with gu we expect these effects to decrease the effective thermal DM annihilation 
cross section much before the theory becomes non-perturbative. Indeed we find, using the 
scan, an upper bound m Xl < 5 TeV. 

Fig. 4 (upper panels) shows the ratio of masses of the next-to-lightest and the lightest 
FGB, 777,423/777,424, as a function of DM mass m xi for radiatively split DM multiplet (left) 
and in the large mass splitting limit (right). In most of the parameter space satisfying the 
Ddm constraint H 23 is much heavier than A 24 so that the effects of higher resonances are 
indeed negligible. This justifies the use of only the lightest FGB when calculating the DM 
density in the scan. 

Fig. 4 (lower panel) shows the relative radiative mass splitting Arri2\j m xi and Am,3i/m xi 
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as a function of the DM mass, m Xl (both splittings are shown in one plot). In most cases 
the mass splitting is below (9(10%) in order for all x components to lie close to the resonant 
condition, as anticipated. 


2. Scalar dark matter 


For scalar DM the interactions with the visible sector are mainly due to the Higgs- 
portal operator, A h{ 4 >' H), in Eq. (20). The interactions due to the exchanges of 
FGBs are subleading except for the resonant annihilation regions — mjpij'l. By 
adjusting the value of A h one can obtain the correct relic abundance for any mass of 
rn ( j >i irrespective of the lightest FGB mass, rriA 24 , see Fig. 3 (bottom panel). In the cal¬ 
culation of the thermal relic abundance we include the following annihilation channels: 
(p\(j) 1 —> bb, cc, t + t~, W + W~, ZZ, hh and it, see Fig. 6. The annihilation cross sections are 

A 2 H m 2 N c (l — 4m 2 /s) 3 ^ 2 


^(0i0i ->• //) = 


cr(0i0i -> VV) = 

where Cw — 1, cz — 1/2 and 

a(cj)\(f)i —>• hh) = 


8tt^1 - Am 2 Js [{m 2 h - s ) 2 + m 2 r 2 ] 

cy \ 2 h \/l — Amy/s{12my — 4 m\s + s 2 ) 

16,rs d 1 “ 4m */ s [( m S - s ) 2 + m l r l] 

A I yi - 4m|/s[(2m| + s) 2 + mjrj] 


(28) 

(29) 


(30) 


32ws G - Am U s [( m l - s ) 2 + m h r l] 

The thermally averaged cross sections and relic abundances are computed following the 
prescription described in Appendix B. The results of the scan are given in Fig. 3 (bottom 
panel). In Fig. 7 we plot the coupling A h necessary to obtain correct DM relic density 
as a function of the DM mass, rri ( p x . As commented in Sec. Ill B the flavon vevs split the 
lightest DM state (j)\ from the heavier ones, such that only (f)y contributes to Ddm (lower 
dashed line). Instead, if the splitting is too small for 0 2j 3 to decay before freeze-out, all 
three components contribute (upper dashed line). In both cases, requiring the Higgs-portal 
coupling A H < \/47r, such that the relic-abundance calculation is well in the perturbative 
regime, limits the DM mass m ( p l < 8 TeV. 

Note that the role of the Higgs portal may be played by other light scalars. I 11 Ref. [20] 
the flavon field of the Abelian horizontal symmetry was used to enhance the DM annihi¬ 
lation cross section. If the flavons are light, they can also modify the phenomenology of 
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Figure 7: The Higgs-DM coupling, Ah, as a function of DM mass that gives the correct relic 
abundance for the Higgs portal scalar DM (red band). The upper (lower) dashed edge corresponds 
to the limit where 02,3 decay much after (before) the thermal freeze-out of 0i. The LUX bound, 
assuming correct relic abundance, is shown as a shaded grey region. 

the fermionic flavored DM, allowing DM annihilation into flavons. I11 this case the DM 
phenomenology of the fermionic flavored DM would be closer to the one of our scalar DM 
model. 

C. Cosmology 

The heavier flavored DM states, both for the fermionic DM, X 2 , 3 , and scalar DM, 02,3, are 
unstable. They decay through the Xi XjQQ' transition when the mass splitting is larger 
than the pion mass, and through the Xi ~ ► Xjll otherwise, cf. Eqs. (17), (18). The SM 
particles in the final state of these decays can have various observable effects in cosmology 
and astrophysics. 

The two relevant sets of parameters are the lifetimes of the two heavy states, r X2 3 , and the 
related mass splittings of the DM multiplet (with respect to the lightest state), Am^i, Am 2 i. 
The lifetimes determine at which cosmological epoch the heavy states decay. The mass split¬ 
tings control the released combined electromagnetic and hadronic energy, E vis ~ Am 2 i, 3 i. 
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They also determine the relic abundances of the heavy states. Generically, near the degen¬ 
erate limit each state contributes roughly a third of the total DM relic abundance, Ddm- 
Close to the resonant condition m x ~ rriA^/2 the Xi, 2,3 relic abundances may differ from 
^dm/ 3, depending on the common DM mass and relative mass splittings. 

For the scalar DM the mass splitting is expected to be large. The 02,3 therefore decay 
before primordial nucleosynthesis. The decays yield negligible entropy release due to the 
small 02,3 abundances. Such scenarios are basically unconstrained by current cosmological 
observations. The same is true for the fermionic DM if additional spurions split the DM 
multiplet at tree level. 

If the fermionic DM multiplet is split solely by radiative corrections, the ^ 2,3 and Xi are 
generically much more degenerate, cf. Fig. 4 (right). The ^ 2,3 states are then potentially long 
lived. For 72,3 ~ (9(10 _1 s — 10 12 s) the decays may affect the primordial generation of light 
nuclear elements [29]. For longer lifetimes, r 2) 3 ~ (9(10 10 s — 10 13 s), the y 2 ,3 decays distort 
the thermalization of the CMB by injecting high-energy photons into the plasma before 
recombination, which is strongly constrained [30, 31]. The X 2,3 states with lifetimes longer 
than 72,3 > 10 10 s are ruled out, if the injected photons carry energy above the thresholds of 
the efficient thermalization processes. Typically this is a fraction of the electron mass. For 
even longer lifetimes, r 2 ,3 > 10 13 s, the y 2 ,3 states decay after recombination. This results 
in photons that free-stream to us and can be searched for in diffuse galactic and extra- 
galactic gamma and X-ray spectra. A combination of measurements excludes scenarios with 
72,3 10 26 s all the way down to Am 2 i, 3 i > 0(10) keV [32], 

In the remainder of this section we consider in more detail the region r X2 3 rv_/ 10 _1 s—10 12 s, 
where the big bang nucleosynthesis (BBN) provides the most stringent constraints [29, 33]. 
The injection of energetic photons or hadrons from X 2,3 decays during or after BBN adds 
an additional non-thermal component to the plasma that can modify the abundances of 
the light elements [34—38] . The bounds differ depending on whether the decays result in 
hadronic or electromagnetic showers in the plasma. The most stringent bounds are for a relic 
that produces mostly hadronic showers. This is because the electromagnetically interacting 
particles such as photons and electrons thermalize very quickly by interacting with the tail of 
the CMB distribution until the universe is 10 6 s old. I 11 our case, the decays y 2 ,3 —i > Xi are 
always kinematically allowed for t X2 3 < 10 12 s. The y 2 ,3 decays thus predominantly produce 
a small number of hadronic jets with a combined released hadronic energy f?h a d — A vis . 
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There are three distinct ranges of lifetimes [39]. For 0.1s < t X 23 < 100 s the dominant 
effect is the inter-conversion between protons and neutrons, which overproduces the 4 He 
abundance. For longer lifetimes, 100 s < r X23 < 10' s, hadro-dissociation is the most efficient 
process and the bounds come from the non-thermal production of Li and D. At late times, 
10 7 s < t X23 < 10 12 s, photo-dissociation caused by direct electromagnetic showers or by 
electromagnetic showers from daughter hadrons can lead to overproduction of 3 He. 

We impose the 4 He, D and 3 He constraints 1 using the results in Ref. [39]. The visible 
energy release in the decays is E vis ~ Am 2 i, 3 i. For 100GeV< Am 2 i, 3 i < lOTeV the 
constraints derived from the three relic mass benchmarks in Ref. [39] are well approximated 
by a power-law scaling with E~^\ The exponents for the three constraints are rp He pa 1/3, 
~ 1/2 and 773 n e ~ 1- For inter-conversion and hadro-dissociation the power-law scaling 
is expected to break down at energies below (9(10) GeV due to the presence of hadronic 
thresholds [39]. We thus do not extrapolate the fit results for 4 He and D for Am 2 i i3 i 
below 10 GeV. We assume that the photo-dissociation effects retain approximate power law 
behavior for E vis large compared to the binding energies of the light nuclei, which is of the 
order of few tens of MeV. In our model for r X2 3 < 10 2 F 2 ) s, the mass splitting, Am 2 i i3 i, is 
always above 10(0.1) GeV. Our approximations are thus always valid for ranges of lifetimes 
for which the 3 He constraints are the most stringent. For the deuterium bound, on the other 
hand, the power-law scaling is expected to fail for part of the parameter space where the 
bound is the most stringent, since Am 21i31 can be as low as a few GeV. We have checked 
using the power-law derived bound that these regions are excluded by several orders of 
magnitude. This gives us confidence to conclude that they remain excluded even with a 
more faithful treatment of hadro-dissociation effects. 

In Fig. 5 we show the distribution of \ 2,3 lifetimes in the viable parameter space of the 
fermionic DM model. The cosmological constraints rule out all points with r 2 3 > 100 s, 
which is the range of lifetimes for which the deuterium bound becomes effective. The points 
with lifetimes r 2j 3 < 100 s, on the other hand, are never excluded by cosmological constraints. 
This is the range of lifetimes where the most strigent bound comes from the 4 He abundance, 

1 The measured 4 He abundance has shifted upwards significantly since the publicaton of Ref. [39]. This 
should weaken the constraints for r 2i 3 < 100 s. The upward shift has no consequences for our conclusions 
since we find that the 4 He constraint from Ref. [39] is already never important in excluding the viable 
parameter space in our models. 
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which, however, is not sufficient to exclude any of our fermionic DM model points. 


D. Direct and indirect dark matter searches 


Both fermionic and scalar flavored DM can produce direct detection signal from DM 
scattering on nuclei. For fermionic DM the scattering is due to t-channel exchanges of 
FGBs. For scalar DM the scattering is dominated by the Higgs exchange in the t-channel, 
while the contribution of FGBs is in general negligible. 

The spin-independent interactions with the nucleons for the fermionic flavored DM are 
described by the following effective Lagrangian [40, 41] 

£dir. = / p (x7mX)(P7» + fn{xi»X)(nYn ). (31) 


The Wilson coefficients f p and f n are the couplings to protons and neutrons, respectively, 


fp = £«; 
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(32) 


A 171 


Qy d are the vectorial couplings of FGBs to quarks, defined in Eq. (25). The spin-independent 
DM-nucleon cross section as measured by the LUX experiment [42] is 
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l + =r I ^-l 
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(33) 


where fi xn is the reduced mass of the (y, n) system. The Xenon atomic and mass numbers 
are denoted by Z and A, respectively. We thus have Z = 54, while A varies between 128 
and 134. With the above relations we calculate the DM-nucleon cross section and compare 
it with the current best limits reported by the LUX experiment [42], The results of the 
scan are shown in Fig. 8. Most of the points lie well below the present LUX bound. This 
is a consequence of the fact that the relic abundance is given by the s-channel resonant 
annihilation, while the direct detection scattering is due to t-channel exchanges of FGBs 
and thus not resonantly enhanced. 

For scalar flavored DM the dominant scattering is through t-channel Higgs-boson ex¬ 
change. This leads to the spin-independent scattering on nucleon N = n,p, [43, 44] 
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Figure 8: The predicted spin-independent cross section for DM scattering on nuclei as a function 
of DM mass for radiatively split fermionic DM (left) and in the large mass-splitting limit (right). 
The LUX bound is the brown shaded region. The color coding for the points is as in Fig. 3. 

The Higgs-nucleon coupling is 

f N ’ h = l + lj2rt N) ’ ( 35 ) 

<7 

where the sum runs over the light quarks, q = u,d, s. f g N ^ are defined by the matrix 
elements of the light-quark scalar currents, m^fq^ = ( N\m g qq\N ). For the s quark we use 
the lattice determination fi N) = 0.043 ±0.011 [45]. The matrix elements for u and d qnarks 
depend strongly on 7rlV-scattering data. A Baryon Chiral Perturbation Theory (ByPT) 
analysis of the 7rlV-scattering data gives — 59(7) MeV [46]. This is in agreement with 
a ByPT fit to world lattice Nf — 2 + 1 QCD data, which gives ovv = 52(3)(8) MeV 
[47]. Including both A(1232) and finite-spacing parametrization in the fit shifts the central 
value to oy-Tv = 44 MeV. To be conservative we use oyjv = (50 ± 15) MeV that leads to 
fi p) = (1.8±0.5)-10 -2 , /j p) = (3.4±1.1)-10 -2 , /i n) = (1.6±0.5)-10 -2 , f d n) = (3.8±1.1)-10" 2 
by using expressions in Refs. [48, 49]. This results in = (29.7 ± 1.3) • 10~ 2 where we 
averaged over Higgs couplings to proton and neutron (the difference is an order of magnitude 
smaller than the quoted error). The resulting bound from the LUX experiment, assuming 
correct relic abundance, is shown in Fig. 7 and constrains m ^ > 150 GeV. 

Finally, we discuss the constraints from indirect DM searches. The flavored DM annihi¬ 
lates to quarks so that the most constraining indirect DM searches are due to the photon 
and antiproton cosmic-ray fluxes. The constraints from the antiproton flux are quite depen- 


23 








X 

O H 


.T 

T 

'b 



100 , 


0 . 01 , 


m A 24[GeV] 


> 

* Excluded by H DM 

* Flavor or direct bounds 

* Cosmological bounds 

* All constraints OK 


100 200 


500 1000 2000 


5000 


XI 

a. 


t 

T 

& 

'b 



0.01 j 

ir\-4 • Nonperturbative 

^ • Excluded by H DM 

• Flavor or direct bounds 
IQ-6 _ • All constraints OK 


100 200 500 1000 2000 5000 


m A 24[GeV] 


Figure 9: The dijet production cross section at 8TeV LHC as a function of the lightest FGB mass 
for radiatively split fermionic DM (left) and in the large mass splitting limit (right). The 95% CL 
limit from Ref. [55] is denoted with a solid orange line. The color coding is the same as in Fig. 4. 

dent on the propagation model. This can lead to almost an order of magnitude difference 
in uncertainty on the value of the excluded annihilation cross section [50]. 

For instance, by using the MED propagation model the antiproton-flux measurement 
by Pamela [51] constrains the DM mass to be m x > 20 GeV if the XX' bb annihilation 
dominates. Similar sensitivity is expected from annihilations to other quarks. The FERMI- 
LAT measurements of the photon flux from dwarf spheroidals bound m x > 100 GeV for 
thermal DM annihilating to quarks [52] (there are also slightly less stringent constraints from 
y-ray emissions from the Large Magellanic Cloud [53] , and from isotropic y-ray background 
[54]). 

E. Searches at the LHC 

The searches for particles beyond the SM at the LHC are sensitive to the lightest new 
states in our models. The searches for dijet resonances impose constraints on the mass of 
the lightest FGB [55], and the searches for vector-like T and B quarks impose constraints 
on the mass of the lightest quark partners [56]. 

The FGBs are narrow resonances that have flavor-conserving as well as flavor-violating 
couplings to the SM quarks, ty, d i} and to the quark partners, w', d[. Since the FGBs are not 
colored they do not directly couple to gluons. At the partonic level the production process 
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is dominated by q % q 3 —> A m —> q^qi. The FGBs would then appear as resonances in the 
dijet invariant-mass spectrum. For the most part, the LHC dijet resonance searches are 
relevant only for the lightest FGB which has, to a very good approximation, flavor-diagonal 
couplings to quarks. In this case, the cross section for pp —> jj is given by [57, 58] 

a(pp ^qq) = [ dy B [ dz ./) (y/re yB ) fj (y/re ~ VB ) (36) 
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where, in the partonic center of mass frame, is the total energy, (3f is the velocity of the 
final-state quarks, z = cos 6 * is the cosine of the polar angle of the outgoing quark w.r.t. the 
direction of the incoming quark, and the couplings Q v , Qa of FGBs to quarks were defined 
in Eq. (25). We have only included the s-channel contribution that dominates on the FGB 
resonance peak. Terms odd in z were dropped in the differential cross section since they 
vanish after integration. The predicted dijet cross sections at the LHC with y/s = 8 TeV 
are shown in Fig. 9, where the 95% CL exclusion from Ref. [55] is denoted with a solid 
orange line. This mostly excludes the points where the lightest FGB has large couplings 
to the quarks. Such points are in fact already mostly excluded either by the perturbativity 
requirement or from flavor constraints. 

The quark partners, w', have an inverted mass hierarchy w.r.t. the SM quarks so that 
in most of our scan points the if is the lightest state. The bound on the if mass depends on 
the if —y bW, tZ , and th branching ratios. The respective partial decay widths are given by 

r(T —>• bW ) = \s UL3 V 33 c dL3 \ 2 — y~ (l — 4v) (l + 2a%) > 

047T m W 
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where Xi = m l /m t > and s*, c* are the sines and cosines of the mixing angles between the SM 
and exotic quarks, while V is a unitary matrix describing the misalignment of the Y u and 
Y d vevs. Their definitions can be found in Ref. [23], where also the relevant Feynman rules 
are given. (We present the relevant Higgs Feynman rules in App. C, correcting an obvious 
typographical error of Ref. [23].) In Eq. (38) we took the limit Xb —> 0 that is justified 
since rn t > mb- We use the above expressions for the t! —* bW, tZ , th branching ratios 
to obtain the 95% confidence-level bound on the t' mass by interpolating between quoted 
observed-limits table in [56]. 

F. Flavor constraints 

The model of gauged-flavor symmetries in Eq. (6) was designed to be compatible with 
new TeV-scale dynamics and at the same time satisfy the tight flavor constraints from FCNC 
observables. The FCNCs induced by the exchange of new states are thus relatively mild. 
The light flavor-violating gauge bosons mediate AF = 2 transitions at the tree-level, while 
the light exotic quarks modify the loop-induced SM process. These modifications are large 
enough that they restrict the parameter space of the model [23]. All the flavor-violating 
parameters in the model are fixed by requiring (Y u ) and (Y d ) to reproduce the observed 
masses and mixings in the quark sector. The size of the induced FCNCs thus depend only 
on a relatively small set of flavor conserving parameters in the model, the flavor symmetric 
masses and couplings. Following the analysis in Ref. [23] we focus on AF = 2 observables 
in the neutral B and K sectors, and on B s —y A s q. 

In our analysis we include the mass differences in the neutral A' 0 , B®, and B d sectors, 
A m,K, Am^, and Arris,,, respectively. We also include the indirect CP violation in the 
kaon sector, e K , and the mixing-induced CP asymmetries S^k s and S^. The corresponding 
mixing amplitude 

2 m M (Mg)* = (M\H? s m=2 \M) , (41) 

where M = K°, B d , B®, controls all of these observables. 

Two NP contributions to Mjff dominate. These are the tree-level exchanges of FGBs 
and the loop-induced SM-like contribution with internal up-type quarks, including exotic 
quarks. For the later contribution we first integrate out at the EW scale, jiw — m w, Hm 
exotic quarks together with the W and the top quark. In this step we ignore the hierarchy of 
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masses between the exotic quarks and top. This is a good approximation for the dominant 
contribution that comes from t' . The theory matches onto the EFT with the SM effective 
weak operators. We perform the Renormalization Group (RG) of the Wilson coefficients 
to the low scale at which the hadronic matrix elements are evaluated on the lattice. For 
the tree-level FGB exchanges the hard scale is given by the corresponding gauge-boson 
masses. We integrate out the FGB at the corresponding hard scale and RG evolve the Wilson 
coefficients to the hadronic scale. The FGB exchanges generate four-fermion operators 
with the Dirac structures that differ from the SM one, namely (,/) "f,j. Pr f 3 ) (./)7 ,J ' Pr ./)) and 
( fj ) ( fj'y^PLfi) ) where i,j are the flavor indices. The RG evolution is implemented 
following Ref. [59] (for further details and the dependence of the numerical relevance with 
the scale see also Ref. [23]). For the non-perturbative inputs, the decay constants and the 
bag parameters, we use the current lattice averages [60]. 

The mass difference in the neutral kaon sector, Ai%, and the CP-violating parameter 
£k are given by 


A M k = 2Re(M(f), 


K r e 


IVe 


£k = 


Im (M- 


k°\ 
12 > 


(42) 


V2AM« P 

with ip £ = (43.51 ± 0.05)° and r £ = 0.923 ± 0.006, which includes long-distance effects in 
both Im My^ [61] and in the decay, i.e. I in T [^° [62]. Our SM expectation for Ex incorporates 
the known Next-to-Next-to-Leading-Order (NNLO) QCD corrections due to the charm [63] 
and charm-top [64] contributions. 

The mass differences in the neutral B sectors are given by 


A Ms q = 2|M^ ? |, with cj = d,s. 


(43) 


The CP violation in the neutral B sector is probed by the time-dependent asymmetries in 
the decays B° —> i[>K s and B° s —> i/jcj) that define the observables 


Sq,K 3 = sin(2f3+ 2(p Bd ) and = sin(2|/J s | + 2(f> Ba ) , (44) 


respectively. In the conventional parametrization of the CKM matrix the SM phases are 
given by |e - * 3 and V t ^ M = — | \e~ 1 ^ 3 . The NP phases are defined through 

the relation M^ 1 = |M^ 9 | e 2l ^ q+ ^ B i' ) . The tree level exchanges of FGBs induce such new 
phases in A F = 2 matrix elements. These are thus constrained both by the S^k s and S ^ 
asymmetries, and by ek in the kaon sector. 
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The rate of B —y X s 7 is also modified by the presence of exotic up-type quarks. These can 
only enhance the B —» X s ^ rate with respect to the SM expectations [9]. The contributions 
of FGBs are loop-suppressed. Even though they may be enhanced by rriy they are negligible 
in models with a seesaw-like mass generation for quarks, like the model we consider [65]. The 
SM prediction for the rate in our analysis includes the known NNLO corrections [ 66 - 68 ]. 

In our numerical scan we mark parameter space points to have passed the flavor con¬ 
straints only if the predictions for all our observables lie within three standard deviations of 
the corresponding experimental values. Whenever theoretical uncertainties are relevant, we 
include them in quadrature with the experimental ones. 

The deviations of the selected FCNC observables from the SM predictions for the four 
benchmark points are shown in Figs. 10 to 12. 


V. BENCHMARKS 

To illustrate the most relevant phenomenology of fermionic flavored DM we select four 
representative benchmark points. The main features of the four benchmarks are summarized 
in Figs. 10-13. The upper panels in the figures show the spectra for the FGBs, A m , the quark 
partners, u ' t , d ' t , and the DM multiplet, Xi- Each FGB is represented by four shaded 3x3 
rasters. The shade of the entries in the rasters is approximately logarithmically proportional 
to the size of the couplings to u R , d R , Ul and d,L, respectively (from left to right). The DM 
relic abundances as functions of the Xi mass are shown in the bottom left panels. The 
lines correspond to our approximate calculation for a radiatively split DM multiplet (red 
solid line) and for a DM multiplet with large mass splittings (black dashed line). The 
open diamonds (circles) correspond to the solutions of the coupled Boltzmann equations 
for the radiative (large) splitting cases which were calculated in MadDM. The approximate 
and MadDM relic-abundance calculations are in very good agreement for this small subset of 
benchmarks. In general, however, a disagreement of up (9(30%) could be expected due to the 
approximations (see App. B for a more detailed discussion). The bottom right panels show 
the pull in selected flavor observables, i.e., the differences between theoretical predictions 
and measurements normalized to the 1-er uncertainties. The uncertainties were obtained by 
adding in quadrature the theory and experimental errors. The four benchmark points are 
also marked in Figs. 3, 4, and 5 with a diamond (1), a five-point star (2), a triangle (3), and 
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M„=980.GeV, A*=l., A„ = 11.8, M d = 700.GeV, A d =0.82, A rf =3.8, # e =0.126, g v = 0.41, g fl =0.237 
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Figure 10: Mass spectrum and flavor decomposition (upper panel), DM relic density as a function 
of the DM mass with all other parameters fixed (lower left panel) and the pattern of effects in 
selected flavor observables (lower right panel) for the fermionic flavored DM benchmark 1. The 
input benchmark-point parameters are listed in the center. See text for details. 

a six-point star (4). 

“Benchmark 1” is an example of fermionic flavored DM, where the DM multiplet is 
light, with mass below 1 TeV. The mass of the lightest state in the DM multiplet is m Xl ~ 
520(540) GeV if it lies just below (above) the LFGB resonance. If the mass splitting between 
Xi, X 2 and % 3 is solely due to radiative corrections, y 2 and y 3 are almost mass degenerate with 
masses roughly 10 GeV above m xi , and y 3 is about 100 MeV heavier than y 2 - The lightest 
quark partner is t' with a mass rn t / ~ 1.3 TeV. The lightest FGB has a mass ~ 1.1 TeV. 
All the remaining NP states are above 7 TeV. This benchmark point demonstrates that even 
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parameter regions with low lying FGBs can be consistent with both the resonance searches 
and the FCNC bounds. The most robust constraints in this parameter region come from 
cosmology (in case of radiatively-split DM masses) and dijet-resonance searches (see Figs. 5 
and 9). Note in particular that for the completely (mass) decoupled fermionic DM scenario, 
in which cosmology bounds are absent, all experimental constraints can be satisfied even for 
DM (and LFGB) masses below few 100 GeV. 

The bottom left panel in Fig. 10 shows the predicted relic abundance for this benchmark, 
if only the DM mass is varied, which also modifies the splitting within the DM triplet. Relic 
abundance consistent with observations is obtained for a mass of DM close to half of the mass 
of the lightest FGB, in which case the annihilation cross section is resonantly enhanced. To 
saturate the observed DM relic density, two solutions for m Xl are obtained, with m Xl either 
above or below the resonant peak. We see that for radiatively split DM masses, where all 
components contribute to the DM relic density, m Xi need to he within 0(3%) of the resonant 
peak for the annihilation to be strong enough. For completely decoupled DM multiplet the 
resonant condition is relaxed and needs to be satisfied to 0(5%). 

In Fig. 10 (upper panel) we show the spectrum for the lower mass solution and radiative 
DM multiplet splitting. We see that the quark partners of the lighter generations are heavier 
than the partners of the third generation quarks. Similarly, the FGBs that couple more 
strongly to the first two generations are typically heavier than the ones that couple preferably 
to the third generation. The couplings of the lightest FGB to the light quarks have the form 
Qt — Q'r — Ql — c* -^ 8 j where the relative corrections to this relation are below the 
percent level. This means that the lightest FGB couples to the light quarks vectorially, 
qu 4 ^ g^’ d ^ a very good approximation. The same is true for the majority of parameter- 
space points passing flavor constraints. 

The largest effects in flavor physics are in the mixing observables, the mass splittings 
A md )S in — B,i and B s — B s systems, respectively, and the mass splitting in the K — K 
mixing, Am k , and the related CP violating parameter ex- The pulls in b —> sy and Bd — Bd 
mixing are due to the fact that the measurements agree with the SM prediction only at 1-cr 
level and the contribution to them from new states is very small. 

“Benchmark 2” is an example of a generic parameter-space region, but towards the upper 
end of the perturbatively allowed region. The DM has a mass m xi ~ 4.5 TeV, while the 
heavier states in the DM multiplet have masses 120 GeV and 150 GeV above m xi (for the 
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Benchmark 2 



M„=980.GeV, A„=l„ A„=0.79, M d = 360.GeV, A rf =0.146, A rf =0.91, g e =1.93, g(/=0.47, g D =3A 



Figure 11: Same as Fig. 10 for benchmark 2. 

case of only radiative mass splitting). The lightest exotic quark is the top partner with mass 
rrit> — 1.7 TeV, while the mass of the lightest FGB is ~ 9.2 TeV. 

For such high DM masses it is barely possible to obtain the correct relic abundance (see 
the lower left panel in Fig. 11). Therefore, the DM mass is finely tuned to be exactly on 
the resonant peak (see the lower left panel in Fig. 11), so m Xi ~ 777 ,^ 24 / 2 . Because of the 
high masses of the NP states the direct searches (direct DM detection, t' searches and dijet 
resonance searches) as well as the indirect flavor constraints are easily avoided, although 
K — K mixing does receive non-negligible contributions. 

Also in this case, the couplings of the lightest FGB to the light quarks have the form 
~ Qft ~ Q'l ~ Qft cx A 8 , so the couplings of the lightest FGB to quarks are to a good 
extent vectorial. 
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Benchmark 3 



M„=1210.GeV, A„=l., A„=1.3, M d = 293.GeV, A rf =0.36, A d =1.24, g e =0.45, g v =0A6, g D = 0.82 



Figure 12: Same as Fig. 10 for benchmark 3. 

“Benchmark 3” is an example of a generic parameter space in which all the couplings of the 
model are well below the pertubativity bounds. The lightest FGB has a mass m ^ 24 — 5 TeV 
while all other FGBs have masses above 100 TeV. The lightest partner quark is t' with mass 
m t i ~ 2 TeV. The DM states have masses m xi ~ 2.4 TeV and m X2 ~ m X3 = 2.5 TeV (for ra¬ 
diative mass splitting). All direct experimental constraints as well indirect flavor bounds are 
easily satisfied in this case. For radiatively DM mass splitting the cosmological constraints 
are the most constraining. In particular, requiring small enough r X2 3 (or equivalently large 
enough Am 2 i, 3 i) typically imposes a lower bound on gu . 

“Benchmark 4” is an example of the case in which the next-to-lightest FGBs have masses 
not too far from the lightest FGB. In the benchmark point the lightest FGB has a mass 
m^24 ~ 8.3 TeV, while the next to lightest FGBs have masses 23 ~ 771^22 ~ 19 TeV, 
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M u = 830.GeV, A„=l„ A„=3.5, M d = 1920.GeV, A rf =0.096, A rf =0.48, # e =1.89, g v = 0.312, g D = 1.93 
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Figure 13: Same as Fig. 10 for benchmark 4. 

m^i ~ 39 TeV. In this case the deviations from the ~ ~ Q d L ~ Qf { oc A 8 relation for the 

lightest FGB coupling to quarks are of 0(10%). Nonetheless, this does not have a significant 
effect on the computation of the DM relic abundance. This is demonstrated in Fig. 13 bottom 
left panel where a comparison is shown between the relic-abundance calculation neglecting 
(dashed lines, labelled “Approx”) and including (full lines, using MadDM) flavor off-diagonal 
lightest FGB couplings and contributions of heavier FGBs. For a more detailed discussion 
of these effects see Appendix B. 

The lightest quark partner is t' with mass m t ’ — 1.4 TeV, and is significantly lighter than 
all FGBs and also DM. The DM states are degenerate to a good approximation, with masses 
m Xl = 4.1 TeV and m X2 ~ m X3 ~ 4.2 TeV. Because of the light t' the flavor constraints are 
nontrivial, and there are visible effects in B k and B s mixing observables. 
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VI. CONCLUSIONS 


We investigated the possibility that DM is in a nontrivial representation of the continuous 
flavor group = SU(3)q x SU(3)jj x SU( 3)d- The two main results are that (i) one 
can have a viable model of DM where DM is stable because it is charged under Z 3 - a 
discrete central subgroup of and color £77(3), and (ii) that the DM spectrum can be 
very non-MFV like, while all the low-energy constraints will appear MFV-like. 

Z£ is exactly conserved in many models of flavor. For instance, it remains unbroken for 
MFV new physics. More generally, Z% remains exact if the flavor group t/p M is broken only 
by the vevs of scalar fields, or condensates, with zero flavor triality. Examples of zero flavor 
triality fields are scalars in bifundamental or in adjoint representations of the flavor SU(3) , s. 
The basic requirement for this set-up is that t/p is a good symmetry in the UV. This is 
achieved, if t/p is fully gauged in the UV, which is the possibility we explored. The DM is 
then stable because it is Z£ odd, while all SM fields are Z£ even. 

We investigated two different types of flavored DM models: (i) models in which the 
leading interaction of the DM with the visible sector is through the flavored gauge bosons 
(FGBs), and (ii) models in which the contributions from the FGB exchanges are subleading. 

As an example of the first type of models we considered a Dirac fermion DM that is 
in a fundamental representation of SU(3)u- The relic abundance is fixed by the resonant 
DM annihilation to SM particles through the s-channel exchange of the lightest FGB. The 
DM is thus required to have a mass of about half of the lightest FGB’s mass. This in turn 
implies that the FGBs cannot be arbitrarily heavy, but at most (9(10 TeV). Such light FGBs 
are possible, if the masses of FGBs are inversely proportional to the corresponding quark 
masses. That is, if the FGBs that couple most strongly to light quarks are also the heaviest. 
To achieve this we adopted the model of Ref. [9] in which the inverse proportionality is 
achieved by introducing a set of quark partners also necessary to cancel gauge anomalies. 
The same quark partners also mix with the SM quarks and lead to the mass hierarchy of 
the SM quark masses. 

The flavor and collider phenomenology of the model is very similar to the case where DM 
is not considered. The fact that the first-generation quark partners are the heaviest and that 
the spectrum is completely split, signals the non-MFV character of the model. However, 
the low-energy consequences are MFV-like (see Appendix A). The flavor constraints are 
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satisfied even with FGBs and the top-quark partner, with masses potentially well below the 
TeV scale. The relevant direct collider searches are the searches for dijet resonances and 
t' searches. They exclude part of the available parameter space. Requiring that there is a 
thermal relic DM introduces new constraints. Because DM is part of a flavor multiplet the 
heavier DM states need to decay before big bang nucleosynthesis. In the case of radiatively- 
split fermionic DM this excludes a large part of the parameter space. The remaining points 
are mostly safe from direct-detection bounds. The fact that the DM mass is related to the 
FGB mass by the requirement of almost resonant annihilation sets both lower and upper 
bounds on the DM mass. Requiring that the theory is perturbative also puts an upper 
bound on the DM mass, m Xl < 5 TeV. On the other hand, requiring that the FGBs satisfy 
flavor and direct constraints and that DM is simultaneously in accordance with cosmological 
constraints leads to a lower bound on the DM mass, m Xl > 500 GeV. Improved bounds on 
dijet resonances at the LHC are expected to further strengthen this constraint (see Fig. 9). 

We have also considered the possibility that the DM multiplet is split due to an extra 
source of flavor breaking. Also in this case, the correct relic abundance requires resonant 
annihilation. The DM mass is in thus still roughly equal to half of the mass of the lightest 
FGB. However, the heavier DM states decay well before big bang nucleosynthesis so that 
a much wider range of DM masses is phenomenologically viable. In our scan this includes 
DM masses as light as 100 GeV (with very small couplings to FGBs) and up to 10 TeV. 

A possible signal of the gauged flavor model with fermionic DM at the LHC are mono-jets, 
where the lightest FGB is produced associated with initial state radiation and decays to X\ 
pairs. The yqs are expected to be non-relativistic in the lightest FGB’s rest frame, as their 
combined mass needs to be close to the FGB mass to fulfill the resonance condition for relic 
DM abundance. In the event that such a signal would eventually emerge, the corresponding 
dijet-resonance signal is generically also expected in the model. A final possibility in the 
case of radiatively split DM mass spectrum is that some of the lightest FGBs decay to slow- 
moving X 2 , 3 - They in turn decay within the detector, leaving (highly) displaced vertices, 
isolated hits in the calorimeter or in the muon chambers. LInfortunately, in most of the 
parameter space X 2,3 are expected to decay well outside the detectors, see Fig. 5, leaving 
mono-jets as the only signal. 

In the second type of models, where FGB exchanges give only subleading contributions, 
the only visible consequence of the flavor dynamics on the DM is that DM is stable. The 
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DM mass and the mass of the lightest FGB are no longer connected. We show this in the 
example of scalar flavored DM, in which the dominant interactions with the visible sector 
are through the Higgs portal operator. In this case the phenomenology of the DM is to a 
very good approximation the same as in the Higgs-portal scalar DM, while the dynamics of 
FGBs and quark partners is unrelated to DM. 

In short, we have shown, using an explicit renormalizable model, that it is possible for 
flavored DM to be a thermal relic. The considered model is not the only choice. One 
could consider DM in other representations of C/f. M . Our analysis can be extended also in 
other ways: for instance, by enlarging the global symmetry as in Ref. [13] and subsequently 
gauging it. For instance, with our held content the global group is x SU( 3) x , where 
X is in the fundamental of SU(3) X . In our work we have identified SU( 3) x with SU(3)u, 
but other choices could be made. Yet another possibility is to gauge only part of C/f. M , for 
instance a U( 2) 3 C C/f, M . Note that for fermionic DM, Z 3 is part of an accidental U( l) x 
acting in the dark sector. The U( l) x can be broken by the dimension-7 operator LH\XXi 
but is exact in our renormalizable model. It can in principle be gauged and lead to additional 
phenomenology. If DM is a scalar, U{ l) x can be broken already at the renormalizable level, 
leaving only Z 3 exact. 
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Appendix A: Minimal flavor violation with gauged flavor symmetries 

In this appendix we verify numerically that the Wilson coefficients in the weak Hamilto¬ 
nian for Bd and B s mixing, 
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Figure 14: The FGB contributions to the V—A current operator in the effective weak Hamiltonian. 
The left panel shows the values of the complex ratio 8C\ d /5C\ S for our scan points, with green 
points satisfying all constraints, magenta points excluded by flavor constraints and grey points by 
perturbativity considerations. The cross denotes the point 5C\ d /5C\ S = {V t * d /V t *) 2 . In this point 
the MFV operator with the smallest number of Yukawa insertions completely dominates. The right 
panel shows the cumulative function Pmfv(^), see Eq. (A4). 

generated from exchanges of flavored gauge bosons, are of the MFV type. A tree-level 
exchange of FGBs generates contributions to B c i mixing through operators 

QV = <?“ = 4*5 14K, (A 2 ) 

and Q\ d , Q^ d , that follow from Q bd , Q bd with Pl <h- Pr exchange (the remaining operators can 
be found in, e.g., Ref. [69]). If the lightest FGB has predominantly left-handed couplings, 
then the C\ d Wilson coefficient is the largest one. If the lightest FGB couples predominantly 
to the right-handed quarks, then C\ d dominates. For comparable left- and right-handed 
couplings all four Wilson coefficients, C\ d 3 , C\ d 3 , are important. The analogous discussion 
applies to the case of FGB contributions to B s mixing obtained with a trivial d —)■ s exchange. 

As discussed in Section III the contributions from the gauged flavor model is expandable 
in terms of the SM Yukawas. The contributions due to FGB exchanges can thus be written 
as 

SCY = Ci(mD; 3 + C2(t/„t/bis(wk»l)i3 + • • • = Ci(vsvy 2 + c 2 y 2 d (V^V a ) 2 + ■■■ , (A3) 

where (yd)ij = diag (yd,y s ,yb), and we set y t = 1 in the second equality. In Eq. (A3) we kept 
only the two terms relevant for the discussion below. The same expansion applies for SC bs 
with the replacement d —)■ s in Eq. (A3). 
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In Fig. 14 (left) we show the ratio 8C\ d / 8C\ S , i.e. the NP contribution to the V — A quark 
current operator due to tree-level FGB exchanges. Note that the ratio 8C\ d /8C\ s can be 
complex. If ci, i.e. the leading MFV (LMFV) term, dominates then 8C\ d / 8C\ S ~ iyu/^ts) 2 ■ 
This is denoted by a cross in Fig. 14 (left) . The addition of the operators with extra 
insertions of y d y\ leads to 8C\ d /8C\ s not being equal to (V td /V t * s ) 2 . We verified that the 
curve for 8C\ d /8C\ S shown in Fig. 14 (left) can be fitted with the form of 8C bd ' bs in Eq. (A3) 
by taking Ci real, c 2 complex, and varying Ci from 0(1) to vanishingly small. The points 
in our scan can be grouped into two sets. For the Erst set of points both c\ and c 2 terms 
are sizeable. For the second set of points the ci term is negligible and c 2 dominates. This is 
shown in Fig. 14 (right), where we plot the cumulative distribution function 

-PmfvM =- -r, -• (A4) 

^ total 

The function Pmfv (ti) can be interpreted as the fraction of points that have the ratio 
\8C\ d / 8C\ S \ effectively dominated by operators with up to y d insertions. That is, the points 
dominated by the ci term contribute to Pmfv(O) (and to -Pmfv(^) with n > 0), while the 
points dominated by the c 2 term contribute to P MFV (2) (and to Pmfv(^) with n > 2). The 
points with both Ci and c 2 start contributing to Pmfv(^) for n somewhere between 0 and 2, 
depending on the relative sizes of Ci and c 2 . Fig. 14 (right) shows that the Ci term dominates 
in a subleading (but nonzero) set of points, that about 10% points have both sizeable c\ and 
c 2 terms, and that c 2 dominates in about 80% of the points. 

The similar analysis can be performed for V + A operators, Q hd and Q\ s . We expand the 
FGB contributions to their respective Wilson coefficients in terms of the SM Yukawas 


8c\ d = hiylyuvlyb)^ + C 2 (y\y u yiyb)iz(y\ydy\y u yiyb)iz + ■ ■ ■ 

= c iy 2 d yl(v; d v tb f + c 2 yW b (v; d v tb ) 2 + • • • , 


(AS) 


and similarly for 8C\ S with the d —>■ s replacement. We show in Fig. 15 (left) the ratio 
C\ d /C\ s for our scan. If the c\ term dominates, then C\ d /C\ s = (m d /in s ) 2 (V t * d /V t *) 2 , which 
is denoted by the cross in Fig. 15 (left). The points for which also the c 2 operator (and 
other operators denoted by ellipses above) is important then lie away from the C\ d /C\ s = 
(m d /m s ) 2 (V t * d /V t *) 2 region. We also define a cumulative function 


Pmfv(^) 


N(\sc^/scj-\ > KiMryFgAgP) 

-^total 


(A6) 
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Figure 15: The FGB contributions to V + A current operator in the effective weak Hamiltonian. 
Left panel shows the complex ratio 8C\ d /8C\ S for our scan points with the same color coding as 
in Fig. 14. The cross denotes the point 8C\ d /5C\ S = (md/m s ) 2 (V t * d /V t * s ) 2 , obtained if the MFV 
operator with the smallest number of Yukawa insertions dominates. The right panel shows the 
cumulative function PmfvMj see Eq. (A6). 

The values for Pmfv(^) are shown in Fig. 15 (right). We see that also in this case the points 
cluster into two groups, with vanishing c± term or with both c\ and <?2 relevant. 

The above analysis demonstrates that the FGB contributions to the Wilson coefficients 
in the effective weak Hamiltonian can be expanded in terms of the SM Yukawas. This is a 
hallmark of (general) MFV. In particular, the expansion in terms of m^s/mb and off-diagonal 
CKM elements can still be performed and is not ruined by the large ratios of scales present 
in the problem such as the very disparate FGB masses. 

Appendix B: Thermal relic computation 

In this appendix we describe the calculation of relic density that was used in the scans in 
the main part of the paper. Several approximations to the coupled Boltzmann equations were 
necessary in order to reduce the evaluation time per benchmark and thus allow adequate 
coverage of the parameter space. We find the approximate solutions to be in agreement 
with the full solutions at the 0(30%) level. The full numerical solution of the Boltzmann 
equations was obtained with MadDM [27] using a UFO model file [70], which was generated 
with the FeynRules package [28]. 

Denoting the DM multiplet by (p , where (p is either a Dirac fermion or a complex scalar, 
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Figure 16: The fraction of benchmarks as a function of the off-diagonal couplings of the heav¬ 
iest and next-to-heaviest DM components to the lightest FGB (A 24 ) normalized by the diagonal 
coupling of the heaviest component. 


the most general set of coupled Boltzmann equations reads [25] 
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(Bl) 


where X denotes a generic SM state. For large mass splittings between the ip components 
it is sufficient to consider the lightest ipi state in Eq. (Bl). The contributions to the DM 
relic density from the heavy ip components are exponentially suppressed by corresponding 
Boltzmann factors and can be neglected within our precision. In contrast, when the mass 
splittings are small the full set of coupled Boltzmann equations in Eq. (Bl) needs to be 
considered. Nevertheless, even in this case several approximations are possible for our model, 
as we explain below. 

First of all, the coannihilation of different ipi components into SM particles, ifiifj —> XX 
(i t - j), can be safely neglected in our model. In benchmarks that survive the experimental 
constraints the off-diagonal couplings of the lightest FGB to ip are much smaller than the 
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diagonal ones, see Fig. 16. Secondly, in the calculation of DM relic density we also neglect the 
flavor-changing DM scattering off the thermal background, tpiX — y <PjX. The (p t X —> ipjX 
scattering can be important if (a((piX —> <PjX)v i a b) / (a(ipiipj — > XX) n lab ) > ~ 

1CF 9 . In this case the off-diagonal couplings of 0(1CT 4 ) relative to the diagonal ones are in 
principle large enough to have 0(1) effects on the relic density, and neglecting tp^X —>■ ipjX 
may not be justified. Therefore, for the benchmarks with (§^ 23 /(g^f )33 > 3 x 10” 4 and 
small mass splittings among tp, we explicitly verified using MadDM that neglecting ip t X — > <PjX 
scattering leads to a change in DM relic density smaller than 0(30%). 

Finally, we are able to neglect the pure DM scattering process in the third line of Eq. (Bl) 
since (cr(<pi(pj <(->■ <Pk<Pe)v iab) "C (a(ipiipj O XX) t’i a b) in our model. The largest contribution 
to this process is from diagonal couplings between the FGB and DM. This process can cou¬ 
ple the evolution of the DM species if ( o(<piipj <B- ipk<pc)v i a b) ~ ^ XX) v^). The 

diagonal FGB couplings to the quarks and to the DM of the same generation are approx¬ 
imately equal. By accounting for the color factors and the multiplicity of channels when 
annihilating into SM fields one concludes that the pure DM scattering is indeed subleading. 

Therefore, for almost mass degenerate ipi it is sufficient to consider a set of uncoupled 
Boltzmann equations 

+ 3 Hn Vi =-^2 (viViVi +* XX ) Dab) (nj. - n^f) . (B2) 


o 

The DM relic abundance is in this case the sum of relic abundances for each of the three 
components obtained from the above set of equations (the heavy (p components, in our 
case ip 2 and <p 3 , decay after their respective freeze-outs and contribute to the (pi DM relic 
abundance). In contrast, for large mass splittings the heavy (p components are irrelevant 
for the calculation of the DM relic abundance. This is then obtained from Eq. (B2) by 
considering only the lightest DM state, in our case <pi. 

We calculate the DM relic abundance from Eq. (B2) using the freeze-out approxima¬ 
tion [26], which gives 


(B3) 

J (xf)y/g^Mpi 

where Mpi = 1.22 x 10 19 GeV, g * is the total number of effectively relativistic degrees of 
freedom at the time of the freeze-out, and 



(B4) 
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The freeze-out temperature (xf — m ip jTf') is obtained by solving 


Xf = In 


0.038 g eS Mpi m vi (av Uh ) 


/th 


A /9* x f 


where the thermally-averaged cross section is 


( cry lab) t h — 


2x2 


7T 


ft <’ ” * , 


(B5) 


(B 6 ) 


with t>i a b — 2\/e(l + t)/(l + 2e) and t — s/{2m lpi ) 2 — 1. The freeze-out approximation 
is accurate to a few percent with respect to the full numerical solution of the Boltzmann 
equation [25]. 

The fermionic flavored DM annihilates through the s-channel exchange of FGBs. In this 
case, the integration over x can be performed analytically and the double integral in Eq. (B4) 
reduces to a single one that can be efficiently evaluated numerically. In particular, 


J(xf) — / 2crui a bErfc (y/xje) de . (B7) 

Jo 

We evaluate the above integral numerically in the parameter scan. 

In the annihilation cross section of the fermionic flavored DM we keep the dominant 
contribution - the s-channel exchange of the lightest FGB, A 24 . The annihilation cross 
section for XiXi UjUj (and similarly for XiXi ~* djdj ) is given by 


• \ \ ' UjUj) 



(B 8 ) 


where the vector and axial-vector couplings to quarks, Ga,v, were defined in Eq. (25), a/s is 
the center-of-mass energy and T ^4 is the total decay width of the lightest FGB. The decay 
rate for A 24 UjUj assuming 24 > 2m Uj is 


r(V 4 -> UjUj) = 


m 


24 


4m 2 . 


47T 


m 


^24 


(9 


u\2 

V 'jj, 24 


1 + 


2 ml 


m 


A 24 


+ (G 


u\2 

A 'jj, 24 


4m 2 


m 


A 24 


(B9) 

The rate for A 24 —> XiXii djdj is obtained after trivial replacements for masses and couplings 
(and dividing by the N c color factor for decays to XiXi)- The total FGB decay rate is 
obtained after summing over all kinematically allowed decay channels. 
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Appendix C: Higgs coupling Feynman rules 


As noted in Sec. IV E, the hff Feynman rules given in Appendix A.l of Ref. [23] contain 
a typo. The corrected Feynman rules are given here. 



where the couplings Cl and Cr are: 


i 
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{C L P L + C R P R ) 


huiUi 

• Cr Cr S UR i C-uli 

hu'iu'i 

• CL Cr X u C UR i Su L i 

huiu[ 

J Cr X u C UR i C UL i 


S UR i Suli 

hu'jUi 

J CR X u S UR i s UL i 


^ Cr X u C UR i C UL i 


(Cl) 


(C2) 
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